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Abstract. 

Between 2000 June and September an exceptional variability phase 
occurred in the peculiar F-type hypergiant p Gas, when the R-brightness 
dimmed by at least a full magnitude. The star recovered from this deep 
minimum by 2001 April. It is the third outburst of p Gas on record in the 
last century. We observe TiO absorption bands in high-resolution optical 
and near-IR spectra obtained with the Utrecht Echelle Spectrograph dur¬ 
ing the summer of 2000. TiO formation in the outer atmosphere occurred 
before the deep brightness minimum. Atmospheric models reveal that the 
effective temperature decreases by at least 3000 K, and the TiO shell is 
driven supersonically with M~5.4xl0“^ Mq yr~^. Strong episodic mass 
loss and TiO have also been observed during the outbursts of 1945-47 
and 1985-86. 

An analysis of the exceptional outburst spectra of 2000-01 is pro¬ 
vided, by comparing with high-resolution optical spectra of the early M- 
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type supergiants // Cep (la) and Betelgeuse (lab). During the outburst, 
central emission appears above the local continuum level in the split Na 
D lines. A prominent optical emission line spectrum appears in variabil¬ 
ity phases of fast wind expansion. The radial velocity curves of Ha, and 
of photospheric metal absorption lines signal a very extended, and ve¬ 
locity stratified dynamic atmosphere. The outburst spectra indicate the 
formation of a low-temperature, optically thick circumstellar gas shell 
of 3x10“^ Mq during 200 d, caused by dynamic instability of the up¬ 
per atmosphere of this pulsating massive supergiant near the Eddington 
luminosity limit. 

We present an equation that correctly predicts the outburst time- 
scale, whereby the shell ejection is driven by the release of hydrogen 
recombination energy. The observations reveal that during the outburst 
the atmospheric hydrogen recombination zone dissipates ~5 % of the 
stellar luminosity into the circumstellar environment, by driving a su¬ 
perwind. We observe that the mass-loss rate during the outburst is of 
the same order of magnitude as has been proposed for the outbursts 
of rj Carinae. These phases of punctuated mass-loss represent the ma¬ 
jor mass-loss mechanism of cool massive hypergiants. The research re¬ 
sults in this paper are described in further detail in Lobel et al. (2003). 
A spectral movie sequenc e of the outburst is available at http://cfa- 


WWW, harvard.edu/^alobe]. 


1. Introduction 

Cool hypergiant stars as p Cas (HD 224014) and HR 8752 are thought to be post- 
red supergiants, rapidly evolving toward the blue supergiant phase (de Jager 
1998). They are rare enigmatic objects, which we are continuously monitoring 
with high spectral resolution since about one decade. Yellow hypergiants (la"*") 
are very important objects to investigate the physical causes for the luminosity 
limit of evolved stars. A spectroscopic study during 1993-95 showed that p Cas is 
a slowly pulsating supergiant with changes in the optical spectrum corresponding 
to variations in Teg of less than 750 K (Teg=7250 K — 6500 K; Lobel et al. 1998). 

During a famous outburst of p Cas in 1945-47 many zero-volt excitation en¬ 
ergy lines appeared, not previously observed in its absorption spectrum. These 
atomic lines, normally observed in M-type supergiants, were strongly blue- 
shifted, signaling the ejection of a cool circumstellar gas shell (Beardsley 1961). 
In the following two decades after this event, several papers were published in 
the astrophysical literature, discussing p Cas’ peculiar spectroscopic changes in 
the years during and after this dramatic outburst (e.g. Popper 1947; Thack¬ 
eray 1948). Between 1945 and 1946 the star rapidly dimmed and the spectrum 
developed optical and near-IR TiO bands. Within a couple of years, the hyper¬ 
giant brightened up by nearly a full magnitude, and a mid G-type spectrum was 
recovered around 1950. More recently, Boyarchuk, Boyarchuk, &: Petrov (1988) 
also reported the detection of TiO bands in the optical spectrum during a more 
moderate outburst of 1985-87. 
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Figure 1. The F-brightness curve of p Cas is compared in the upper 
panel to the radial velocity curve, observed over the past 8.5 years. Ob¬ 
servation dates of echelle spectra are marked with short vertical lines in 
the lower panel. The radial velocity curve of Fe i A5572 shows a strong 
increase of the photospheric pulsation amplitude before the outburst of 
fall 2000 (JD 2451800—JD 2451900), when TiO bands develop {marked 
TiO). 

This paper discusses the spectroscopic changes we monitored during the 
2000-01 outburst. We analyze prominent TiO bands newly detected in the 
spectra, before and during the outburst minimum. A brief discussion of the 
stellar outburst physics is also provided. 

2. Radial Velocity and Brightness Cnrves 

The upper panel of Figure 1 shows photo-electric observations of p Cas in the 
F-band {black dots) by Percy, Kolin, &: Henry (2000) over the past decade, sup¬ 
plemented with visual magnitude estimates from AAV SO, and from AFOEV 
(French Association of Variable Star Observers) during the outburst of late 2000 
{cyan dots). The brightness curve shows semi-regular variability, with the deep 
outburst minimum of F~5™.3 in 2000 September-November, preceded by a 
conspicuously bright visual maximum (F~4™) in 2000 March. 

The lower panel shows the radial velocity curve which has been monitored 
from the unblended Fe i A5572 absorption line {black dots). The radial velocity 
curve, determined from a linear interpolation of the temporal Fe i line profile 
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Figure 2. Near-IR TiO absorption bands are observed during the 
outburst of p Cas in 2000 July. The weak TiO bands around 7070 A 
have distinct shapes, also observed in Betelgeuse. The synthetic spec¬ 
trum of p Cas is computed with only TiO lines, while that of Betelgeuse 
also includes molecular and atomic lines. 


changes, is compared with the F-magnitude curve in the upper panel {blue dotted 
line). We observe that the star becomes brightest for variability phases when 
the atmosphere rapidly expands. J^ad decreases by ~20 kms“^ in less than 200 
d during the outburst event. The short black vertical lines mark a total of 78 
echelle spectra observed with high-resolution spectrographs over the past 8.5 
years. 

The spectra have been obtained from our long-term monitoring campaigns 
with four telescopes in the northern hemisphere; the Utrecht Echelle Spec¬ 
trograph of the William Herschel Telescope (La Palma, Canary Islands), the 
Sofin spectrograph of the Nordic Optical Telescope (La Palma, Canary Islands), 
the Ritter Observatory telescope (OH, USA), and the Zeiss-1000 telescope of 
the Special Astrophysical Observatory of the Russian Academy of Science (Ze- 
lenchuk, Russia). The spectra marked with a colored long vertical line {labeled 
1 to 5) in Figure 1 are discussed below. 

3. Formation of Titanium-oxide Bands during the Outburst 

We observe the formation of titanium-oxide (TiO) absorption 7 -system bands 
in the (WHT-UES) outburst spectra of p Cas in July 2000 {red solid line in 
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Figure 3. High-resolution spectra of p Cas around Fe i A6358 show 
that the photospheric spectrum strongly blue-shifts during the out¬ 
burst. Synthetic spectrum calculations indicate that Tefj decreases to 
~4250 K, and the spectrum becomes comparable to the early M-type 
supergiants p Cep and Betelgeuse, with reff=3500 K. A prominent 
emission line spectrum is observed in p Cas during phases of high 
Teff=7250 K, when the atmosphere rapidly accelerates outward. 


Fig. 2). The calculated spectrum with only TiO lines is shown in Figure 2 with 
the red dotted line. A best fit is obtained for a Kurucz model atmosphere of 
Tefj=3750 K and logg'=0. The characteristic TiO bands around 7070 A {marked 
with vertical lines) are also observed in the supergiant Betelgeuse (M2 lab) {black 
solid line), computed with TpflF=3500 K and log 0 =—0.5 {black dotted line) (see 
Lobel & Dupree 2000). 

The spectrum synthesis typically contains ~1500 TiO lines per A from 
the P, Q, and R-branches, for five ^Ti^^O isotopomers (Z=46-50), with Earth 
abundance fractions (Schwenke 1998). The TiO bands in p Cas signal the for¬ 
mation of a cool, optically thick, circumstellar gas shell with r<4000 K, caused 
by the supersonic expansion of the outer atmosphere during the outburst. We 
observe an expansion velocity of Uexp=35±2 kms“^ for these TiO bands, about 
15-20 kms“^ faster than the atmospheric expansion velocity we determine from 
the atomic photospheric absorption lines. 
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Figure 4. The Ha absorption core becomes very weak during the 
outburst. The line core does not blue-shift, and develops strong emis¬ 
sion line wings. The line becomes hlled in by recombination emission 
with the cooling of the Ha-envelope during the outburst. Strong emis¬ 
sion is also observed in the blue wing of Ha for the pre-outburst cycle, 
during the atmospheric collapse that precedes the large H-brightness 
maximum before outburst. The Ha line core of p Cas is much weaker 
than in Betelgeuse and p Cep, where it is excited by a permanent 
chromosphere. 


4. Modeling the Outburst Spectrum 

In Figure 3 a best fit to the atomic spectrum, observed during the outburst 
of p Cas {solid red line labeled 4), is computed with reff=4250 K and log5=0 
{red dotted line). The graph is centered in the heliocentric velocity scale around 
the Fe i and V i blend at A6358. The vertical dotted line is drawn at the stel¬ 
lar rest velocity of —47 kms“^ (Lobel 1997). During the outburst, the entire 
photospheric spectrum Doppler shifts toward the shorter wavelengths, and de¬ 
velops many atomic absorption features also observed in Betelgeuse and p Cep 
(compare with lower black lines). 

We also observe many peculiar metal emission lines appearing above the 
stellar continuum level during two phases with very fast atmospheric expansion 
in 1993 December {blue vertical line labeled 1 in Figure 1) and 1998 October 
{labeled 2). The accelerated stellar wind collides with circumstellar or interstellar 
material, which excites the permitted emission line spectrum. Further analyses 
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Figure 5. During the outburst central emission appears above the 
local continuum level in the Na -D lines. The broad absorption por¬ 
tions of the lines are strongly intensity saturated, and do not show 
Doppler shifts with time. The maximum of the central emission line is 
red-shifted with respect to the stellar rest velocity, and a blue-shifted 
narrow feature is observed in both lines of the doublet. 


of the formation of the peculiar emission line spectrum in p Cas are presented 
in Lobel (1997 Chapter 4; 2001a). For these variability phases we compute that 
Teff=7250 K {upper blue dotted eurve in Figure 3). The spectrum modeling 
therefore reveals a decrease in of at least 3000 K during the outburst of p 
Cas in late 2000. 


5. Mass-loss Rate During the Outburst 

We compute from the synthetic spectrum fits to the TiO bands a mass-loss rate 
of M~5.4 X 10“^ Mq during the outburst of p Cas. We assume spherical 
geometry; M=47ri2^/9Uexp, with i2*=400±100 Rq, and p~10“^° grcm“^ the 
minimum density in the TiO line formation region. The spectrum synthesis 
reveals that Teff returned to ~5750 K within 100 d after the deep outburst min¬ 
imum {magenta dotted line in Figure 5), corresponding to a brightness increase 
of ~0™.5 {magenta vertieal line labeled 5 in Figure 1). 

During the outburst, the Ha line profile develops prominent emission wings, 
signaling a large M, around a central absorption core that becomes unusually 
weak {red line in Figure 4)- The comparison with the Ha profiles of Betelgeuse 
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and /i Cep in Figure 4 {lower black lines) reveals that the line core during the 
outburst is not excited by a stellar chromosphere, as is the case for these M-type 
supergiants. The spectrum computed for Betelgeuse without the model for the 
chromosphere {black dotted line) shows that the Ha transition is not sufficiently 
excited to match the depth of the observed absorption line core. The Ha line 
profile of p Cas is filled in by recombination line emission during the outburst 
cooling of the upper atmosphere. The line thereby assumes a self-absorbed 
shape, indicating a very extended Ha envelope which is optically thick. In 1999 
May we observe {green vertical line labeled 3 in Figure 1) a very prominent blue 
emission wing in Ha, while the photospheric spectrum strongly shifts to the 
longer wavelengths {green line in Figure 4)- It indicates an exceptional (pre¬ 
outburst) variability phase with very strong downflow in the upper and lower 
atmosphere, which precedes the brightness maximum of 2000 March, leading up 
the actual outburst event. 

6. Na D Emission During the Outburst 

During the outburst we observe prominent emission lines in the deep and broad 
absorption cores of the split Na D doublet lines. The central emission maxima 
appear above the local continuum level in Figure 5 {red lines labeled 4)- This 
results from the strong decrease of the optical continuum flux due to the decrease 
of Teff by at least 3000 K. We observe that the TiO bands develop before the deep 
brightness minimum is reached, when the photosphere assumes the maximum 
expansion velocity in Figure 1 {vertical red line labeled 4)- The H-brightness 
minimum results from a strong decrease of the entire atmospheric temperature 
structure. Unlike the photospheric absorption lines, the Na D lines are intensity 
saturated and do not reveal Doppler shifts with time {compare with 1 and 3 
in Figure 5). The far violet extended Na D line wings form in an optically 
thick wind during the fast atmospheric expansion. These extended line wings 
were very strong in 1993 December {blue line labeled 1), with maximum wind 
velocities up to —130 kms“^. 

7. Dynamic Spectra of 1993 2002 

Figure 6 shows the dynamic spectra of Ha, Fe i A5572, and the split Fe i A5506 
line, observed between 1993 November and 2002 February. Red color indicates 
the stellar continuum level, and blue color the flux depression in the absorption 
lines. The white spots in Ha are emission above the stellar continuum level. 
The line profiles are linearly interpolated between consecutive observation nights, 
marked with the left-hand tickmarks. The radial velocity curves at half intensity 
minimum of Ha, and of the Fe i lines {white dashed lines) reveal a strongly 
velocity stratified dynamic atmosphere. We observe that Ha varies over a much 
longer period of time compared to the photospheric Fe i lines. It signals that 
Ha forms on average higher, over a more extended portion of the pulsating 
atmosphere. Notice the strong blueshift of the Fe i lines during the outburst of 
mid 2000. The outburst is preceded by very strong line emission in the short 
wavelength wing of Ha, while the absorption core extends longward, and the 
photospheric Fe i lines strongly red-shift. We observed that a strong collapse of 
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Figure 6. Dynamic spectra of Ha, Fe i A5572, and the split Fe i 
A5506 line. The white spots in Ha are emission above the stellar con¬ 
tinuum level. The line profiles are linearly interpolated between consec¬ 
utive observation nights, marked with the left-hand tickmarks. Time 
runs upward, indicated for each new calendar year with the right-hand 
numbers. The dashed white lines trace the radial velocity determined 
from the lines. The curves reveal a strongly velocity-stratified dynamic 
atmosphere. Notice the strong blue-shift of the Fe i lines during the 
outburst (mid 2000). The outburst is preceded by very strong emission 
in the short-wavelength wing of Ha, while the absorption core extends 
longward, and the photospheric Fe i lines strongly red-shift. A strong 
collapse of the entire atmosphere precedes the outburst. 
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the entire atmosphere precedes the outburst event during the pre-outburst cycle 
of 1999. 


Notice further the development of even stronger emission in the short- 
wavelength wing of Ha, and very weak Fe i lines during the more recent bright¬ 
ness maximum of 2002 December. High-resolution observations of 2002 June 
and July reveal a remarkable split Ha absorption core. This type of line shape 
was not observed in the past 8.5 years. Only Beardsley (1961) mentions the ob¬ 
servation of “core emission” in Ha during the famous outburst of 1945-47. Our 
Ha observations may therefore signal that a new, presumably stronger, outburst 
is imminent in p Cas. 

A spectral movie of the p Cas outburst is electronically available at 
http://cfa-www.harvard.edu/~alobel through animated GIF and AVI for¬ 
mat files. A high-resolution atlas of the optical spectrum of p Cas will also be 
posted at this URL address. 


8. Hydrogen Recombination Driving Mechanism 

We compute that the hypergiant atmosphere becomes dynamically unstable be¬ 
cause the first generalized adiabatic index Fi=((ilnP/dln/o)ad) volume averaged 
over the atmosphere, assumes values below the stability value of 4/3 (Lobel 
2001b). <Fi> <4/3 when increases to >7250 K during the strong atmo¬ 
spheric collapse before the outburst. The outburst event, directly following this 
strong collapse, is driven by the release of ionization energy due to the recombi¬ 
nation of hydrogen with the cooling of the entire atmosphere by Areff~3500 K. 
For an adiabatic spherical expansion we calculate that the outburst time-scale 
is correctly predicted with (Lobel et al. 2003): 


R* 1 

"" 3ViVfcTo IW 





TeS 



where 



27 + 2(1 + 1 ^, 


21 + 2p^ 


( 1 ) 


Fs is the third adiabatic index, evaluated for a mean ionization fraction 
of 1/2, in the partial hydrogen ionization zone, where the atmosphere is most 
expandable (or compressible). For R*=400 Rq, Tef[=7250 K, and ro~8000 K in 
the partial hydrogen ionization zone with 7=13.6 eV, we compute fburst=221 d, 
and an expansion velocity of 36 kms“^, in good agreement with the observed 
values. Hence, equation (1) predicts that during the outburst the atmosphere 
expands over ~2.5 7?* with the complete recombination of hydrogen. 

We also estimate that the time-scale over which the total hydrogen recom¬ 
bination energy, released into the circumstellar environment by the outburst, 
is dissipated by the stellar radiation field during a period following the fast 
envelope expansion of 

INpds 


^brightness—min 
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where a is the Stefan-Boltzmann constant. For an atmospheric expansion of dg 
= 869 Rq during outburst, and a mean envelope density of p = 10“^° grcm“^, 
we compute with equation (2) that ibrightness-min—86 d for Teff=3500 K. This 
time-scale corresponds very well to the period over which the deep brightness 
minimum of ~100 d was observed in late 2000. The total energy released into 
the circumstellar environment Flsheii = 4:Tt R'^dgpN I with complete hydrogen 
recombination is ~6.1 xlO^^ erg. This amount of energy yields a luminosity of 
Lgheii = 1.8 X 10^ Lq, when radiated during 100 d. For a total stellar luminosity 
of L*= 3.5 X 10^ Lq during quiescent variability phases, Lgheii is ^5% of L*. 

In conclusion, we find that during the outburst the partial hydrogen re¬ 
combination zone of the hypergiant atmosphere dissipates several percent of the 
stellar luminosity into the stellar environment, thereby driving an expanding 
wind to highly supersonic outflow velocities during a global cooling of the entire 
atmosphere by more than 3000 K. 


9. Conclusions 

We observe the formation of TiO bands in the spectrum of p Cas before the 
deep brightness minimum in the outburst of 2000-01. A supersonic expansion 
velocity is observed for the new TiO bands, while V rapidly dims by more 
than a full magnitude. The TiO shell expands faster than the photosphere, 
observed from atomic metal absorption lines. A large oscillation cycle and a 
very bright light maximum, with above 7250 K, precede the outburst event. 
We compute that M~5.4 x 10“^ MQyr~^ during the outburst, whereby T^g 
decreases from ~7250 to <3750 K. T^g returns to 5750 K within 100 d after 
the deep outburst minimum. Since recurrent outbursts occur about every half 
century in p Cas, these outburst phases of punctuated mass-loss are the major 
mass-loss mechanism of this massive cool hypergiant. A prominent emission line 
spectrum is observed above the stellar continuum level in p Cas during phases 
of fast wind expansion. We compute that the outburst is driven by the release 
of hydrogen recombination energy with the cooling of the entire atmosphere. 
Based on observations of the Ha line profile evolution in the summer of 2002 we 
expect that a new outburst of p Cas is imminent. 
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